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Abstract. A one-dimensional system of masses with nearest-neighbor interactions and 
periodic boundary conditions is used to study mode decay and ergodicity in nonlinear, dis- 
ordered systems. The system is given an initial periodic displacement, and the total system 
energy within a specific frequency channel is measured as a function of time. Results indicate 
that the rate of mode decay at early times increases when impurities are added. However, 
for long times the rate of mode decay decreases with increasing impurity mass and impurity 
concentration. This behavior at long times can be explained by Anderson localization effects 
and the nonergodic response of the system. 
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1 Introduction 

The transition from quasi-periodic to ergodic behavior in nonUncar systems has been 
an active field of research since the seminal work of Fermi, Pasta, and Ulam (FPU) [Q . 
The behavior of the FPU system has been discussed in great detail |^ , and given the 
degree of nonlinearity and initial energy density, one can estimate the time required 
before the system becomes ergodic. However, little is known about the effects that 
impurities have upon the transition to ergodic behavior in nonlinear systems. Since 
mode coupling in nonlinear systems occurs through interactions, one expects that the 
presence of impurities will hasten the transition to ergodic behavior. This expectation 
is, however, in general not correct due to Anderson localization effects 

Reported here are the results from a numerical experiment using a one-dimensional 
system of masses with nonlinear nearest neighbor forces and periodic boundary condi- 
tions. The displacement of each mass is sampled over a finite interval of time and the 
energy within all frequencies is calculated for each mass. The total energy in a single 
mode is conserved in harmonic systems both with and without impurities. Therefore, 
for the nonlinear systems, the effects of the impurities upon the mode decay can be 
compared directly to the systems without impurities. 
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2 Numerical Experiment 

The system used in this experiment, which corresponds to the FPU /3-model, is com- 
posed of N masses with periodic boundary conditions, and with unit equilibrium 
spacing. The masses undergo a displacement u{t), and the i-th mass interacts 
through nearest neighbor forces: 

mjitj = -K [(Ui+i - Wj) - (ui - - 13 - u^)^ - (u,, - (1) 

For this experiment, the coupling constants n — P — I. A pure system is composed 
of masses m° = 1, and disorder is achieved by randomly changing a number Nj of 
the masses to a second value. The time integration is performed using a 4-th order 
predictor-corrector algorithm. Further details will be given elsewhere. 

The initial condition is a unit amplitude, zero- velocity sinusoidal displacement with 
wavelength A = 32. Time is scaled by the harmonic frequency uj" = 27r/A. At certain 
intervals, the time-dependent displacement of each mass is transformed to frequency 
space, and these data Uk{xi) are used to calculate the modal energy — E{uo°) and 
mass energy Ei = E{xi) by summing over masses and frequencies, respectively. 



3 Results 

The systems consist of = 256 masses, which is a size consistent with systems 
used elsewhere Time is expressed as the dimensionless quantity uj°t/2'K which 
is equivalent to the number of harmonic cycles, and the modal energy is scaled to 
E{ijj°t/2Ti — 32). The error bars shown in the figures represent the estimated standard 
deviation in the mean. 

The initial systems investigated contained impurity masses M/ > m° . Comparisons 
of the effects upon modal decay as a function of either impurity mass Mi or number of 
impurities Nj are shown in Figs. |l| (a) and (b); the pure nonlinear system is denoted by 
open circles. As expected, the addition of impurities increases the initial rate at which 
the harmonic mode uj° decays. However, for all the combinations of impurity mass 
and concentration shown, the long time decay is slower than for the system without 
impurities. Further, increases in impurity mass and impurity concentration retard the 
long time decay. 

The long time behavior of the systems containing heavy impurities can be explained 
by the nonergodic behavior within the system. Systems containing impurities undergo 
an Anderson transition. The resulting excited modes are localized, concentrating 
energy near the impurities. Since the impurities are heavier than m°, their oscillation 
amplitudes are smaller, resulting in a smaller nonlinear contribution to the energy at 
the impurity. 

As a demonstration of the localization of energy at the impurities, the mass energy 
Ei is calculated at each impurity. The average energy at an impurity {Ei)i is compared 
to the average energy at all the masses {Ei)^- The ratio a — {Ei) i / {Ei) ^ , which is 
interpreted as a measure of the extent to which ergodicity has been achieved, is shown 
in Figs. 0(c) and (d). 
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Fig. 1 Effect of impurity mass Mj and number of impurities Ni on the mode energy E^j 
and relative energy a at an impurity. In (a) and (c) the number of impurities Ni is 4. In (b) 
and (d) the impurity mass Mi is 7. The open circles are for the nonlinear system without 
impurities. Error bars represent the estimated standard deviation of the mean. 

Figures 0(c) and (d) show that the relative energy at an impurity increases with 
increasing impurity mass, and is somewhat insensitive to changes in the impurity con- 
centration, respectively. Therefore, the total energy located at impurities is propor- 
tional to the mass and the number of impurities. This explains why the modal decay 
rate decreases as either the impurity mass or the impurity concentration increases. 

Based upon the arguments given above, impurities with masses that are lighter 
than m° should have a different effect upon mode decay. Since the impurities are 
lighter, the localized energy will create large oscillations, which should enhance non- 
linear interactions and, hence, enhance the rate of mode decay. The results of that 
experiment for impurity mass Mj = 0.1 are shown in Fig. ^ For low concentrations, 
the mode decay shown in Fig. |^(a) is faster than for the pure system, as expected. 
However, at sufficiently high concentration, the mode decays slower than the pure 
system when the average spacing between impurities is one half the initial wavelength. 

The relative energy a located at a light impurity shown in Fig. ||(b) is less than 
unity. The light impurities apparently expel energy due to the large oscillation am- 
plitudes, hastening mode decay. Above a critical concentration, it is conjectured that 
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Fig. 2 The modal energy E^, and the relative energy at the impurities a for impurities of 
mass Mi = 0.1. 



the light impurities interact in such a way as to prolong the mode decay. 



4 Conclusion 



For the simulation times studied here, the long time mode decay in disordered non- 
linear systems seems to be controlled by disorder and Anderson localization effects. 
For either heavy or light impurities the energy becomes localized, and the assumption 
of ergodic behavior is not valid. Further, the response of the system depends upon 
whether the impurities are heavier or lighter than the pure system. Heavier impurities 
consume energy and, because of the smaller oscillation amplitudes, release their energy 
through nonlinear interactions very slowly. Lighter impurities appear to expel energy 
through large oscillation amplitudes, hastening nonlinear interactions and mode decay 
at low concentrations, but behave in a manner similar to heavy impurities at high 
concentrations, possibly due to Anderson localization effects. 
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